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SUMMARY 


A general theory of aeroacoustics of homentropic fluid media 
is presented. The definition of sound (unsteady compressible flow) 
and the concept of Bernoulli enthalpy are the fundamental building 
blocks of the theory. It is proved that a Coriolis acceleration 
is necessary to obtain a transfer of energy from the vortical to 
the acoustic mode. The Liepmann pendulum is used to illustrate the 
physical mechanism of noise production by a fluid flow. The basic 
theory is compared with several known aeroacoustic theories. With 
the Lighthill hypothesis, the theory is formally equivalent to the 
Lighthill or Ribner theory in the low Mach number limit. The 
present theory is complete in that all of the basic aeroacoustic 
interact ions can be investigated, in particular, the excitation of a 
fluid flow by sound. 

The interaction of sound with vortex flows is studied in de- 
tail. A formula for vortex core scattering is derived with a plane 
wave analysis. A general theory of acoustic interaction with mult- 
iple weakly interacting line vortices is presented. The theory is 
applied to estimate the Interference noise between an engine and an 
aircraft vortex wake. For the DC-9 In a standard take-off config- 
uration the interference noise is estimated to be 3 or 4 DB at the 
engine peak Strouhal number, for a ground plane observer. 

The noise produced by a corotating vortex pair is calculated 
with the present theory. For low Mach number (M < 0.1) the sound 
power varies as M? in accordance with known results for compact 
2-D flows. For 0.1 < M < 0.3 j the power radiated Is 10 to 15 DB 
less than the compact m 7 law, even though the basic source is 
still quadrupole. For M> 0.3, higher order multipoles contribute 
to the sound power. 

An estimate of jet Impingement noise is given. It is shown 
that the acoustic overpressures are a result of acceleration of tur- 
bulent eddies through the curved flow near the Impingement point. The 
conventional shear noise term Is amplified by a factor of 3 that 
translates into a factor of 9 ( ~ 10 DB) In the sound power. Our 

estimates are in agreement with the order of magnitude of recent 


measurements. The corotating vortex model or its extension to ring 
vortices is suggested as a dynamic model of the jet impingement 
problem. 

The phenomenon of fluid flow excitation by sound is investi- 
gated. When plane waves impinge on the corotating vortex pair, 
resonant excitation or attenuation is possible. The coupling is a 
maximum when the acoustic excitation frequency is approximately 
twice the rotational frequency of the pair. The wave length can be 
quite large compared to the separation of the vortices. Based on 
our simple model, a qualitative explanation of various experimental 
observations with excited jets is offered. 


I. INTRODUCTION 

In Ref. 1, we presented a comprehensive theory of aero- 
acoustics that departs in several fundamental ways from the more 
familiar theories of Lighthill (Ref. 2), Ribner (Ref. 3) and Lilley 
(Ref. 4). Since the original paper was presented, numerous problems 
have been solved (see Ref. 5 and work reported herein). Based on 
these results and discussions with various people involved in aero- 
acoustics research, we have been able to modify and polish the ori- 
ginal concepts. The basic theory and its application to several 
problems of technological importance is the subject of this report. 

In Section II, we give a detailed derivation of the theory for 
the case of constant entropy (homentropic ) flows. The basic con- 
cepts are most easily presented for this case and the resulting 
theory is directly applicable to a variety of important problems. 
Three basic questions of aeroacoustics are posed and discussed in 
the light of the present work and previous theories: 

1. How is sound produced by a "primary" fluid flow? 

2. How is sound processed by a flow? 

3. How does sound affect a primary flow? 

In Section III, we consider in detail the question of how sound 
is processed by vortex flows. First, we investigate the importance 
of vortex core structure by considering the scattering of plane waves 
from a single vortex. A general formula for the core scattering is 
derived. a Lagrangian approach is used to calculate the scattering 
of an arbitrary sound field from a set of discrete weakly- interact ing 
line vortices. The special case of a line acoustic source is cal- 
culated in detail and used to estimate the scattering of engine noise 
from an aircraft vortex wake. Specific application is made to the 
DC-9 where it is estimated that acoustic overpressures of 3 or 4 DB 
can result from the interference between the engine and wake. 

In Section IV, we consider the production of sound by vortex 
flows. First, we consider the model problem of a corotating vortex 
pair and calculate the effects of noncompactness on the sound 
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radiated. Next, we consider the problem of jet impingement noise. 
The shear noise source term is calculated for a turbulent stagna- 
tion point flow. It is estimated that the convective acceleration 
of turbulent eddies In the curved stagnation point flow can pro- 
duce some 10 DB more noise than a free jet. It Is suggested that 
the spinning vortex pair Cor vortex rings') could be used to model 
the impingement problem. 

In Section V, we consider the difficult question of how sound 
affects a fluid flow. The interaction of sound with the spinning 
vortex pair is used to illustrate the basic phenomena. It is shown 
that the vortex motion can be excited or deexcited by sound at the 
acoustic frequency of the vortex (twice the rotational frequency) . 
The results are in qualitative agreement with experimental results 
on excited jets. 
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NOMENCLATURE 
isentropic sound speed 

Lagrange coordinate of line vortex, see Eq. (3.43) 
see Eq. (3-64) 
see Eq. (3.65) 

amplitude of multipoles in vortex noise, see Eqs. 

(4.23) and (4.25) 

decibel 

kinetic energy of primary flow, see Eq. (2.31) 

source function for corotating vortex pair, see 
Eq. (4.11) and Fig. 4.2 

enthalpy 

h + v 2 /2 

Bernoulli Enthalpy, see Eq. (‘2.19) 

Hankel function of first kind 
Heaviside step function 
Cartesian unit vectors, see Fig. 3.1 
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acoustic intensity 

Bessel function of the first kind 

wave number 3 2 tt/X 

Mach number 

pressure 

sound power, see Eq. (*1.29) 

length scale of solid body, exponential and Betz 
cores, see Eq. (3.38) 

unit vector 

see Pig. 4.1 

plane polar coordinates 

normalized radius vector, see Eq. (4.8) 

scattering factor, see Eq. (3.27) 

core scattering factor, see Eq. (3.24) 

see Eqs . (4.18) and (4.19) 

time 

nonacoustic velocity 

mean and turbulent velocity component in a turbulent 
flow 

velocity vector 

acoustic particle velocity 

velocity distribution in vortex, see Eq. (3.1) 
position vector 
see Eq. (5. 7) 

Bessel function of the second kind 


a see Eq. (4.8) 

Y r/2ir modified vortex strength (also ratio of specific 

heats ) 
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total vortex strength, see Eq. (3.25) 
Christoff el symbol, see Eq. (4.38) 

Delta function 

perturbation enthalpy, pressure, density 

see Eq. (5.9) 

kr Q , see Eq. (5.6) 

wave length of sound, 2Tr/k 

vortex acoustic length scale, see Eq. (3.31) 

kinematic viscosity 

density 

dimensionless time 
acoustic potential 
vortex potential 
radian frequency 
curl v , vorticity 
streamf unction 

angular velocity of corotating vortex pair 
vortex core vorticity distribution 


Special Notation 

curl curl operation 

div divergence operation 


D 

Dt 

D 

Dt 

grad 


<q> 


substantive derivative + v«grad 

8 -*■ , 

+ u-grad 

gradient operation 

time average of any quantity q 
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absolute value of vector q 
see Eq. (3-63) 

real part of complex quantity 
imaginary part of complex quantity 


Laplace operator 



II. ACOUSTIC THEORY OP HOMENTROPIC PLOWS 
A. ASSUMPTIONS AND BASIC EQUATIONS 

We consider the class of aeroacoustic problems for which we 
can justifiably ignore real fluid properties that lead to signif- 
icant entropy variations; e.g., the generation of heat by viscous 
dissipation of mechanical turbulent energy and the production of 
entropy by Internal conduction of heat or by nonadiabatic processes 
at a boundary. Clearly, our assumption of homentropic flow imposes 
an upper bound on the Mach number and temperature gradients in the 
flow. Generally speaking we will be concerned with low to moderate 
subsonic Mach numbers and the interaction between two of the funda- 
mental modes of energy transport in a fluid (Ref. 6); 

1 ) The acoustic mode 

2 ) The vortical mode 

Our resulting theory Is directly applicable to the study of sound 
interaction with and production of sound by strong vortical flows. 
The specific application to the scattering of sound by aircraft 
vortex wakes and the production of sound by an impinging jet is 
considered in Sections III and IV. The extension of the basic 
theory to high Mach number nonadiabatic flows will be considered 
elsewhere . 

With the assumption of uniform entropy, the enthalpy, pressure 
and density variations are proportional; i.e.. 



where p is the density and a is the local isentropic speed of 
sound. For a perfect gas with constant specific heats, we have, 
in particular 

a 2 = (y-l)h = (Perfect Gas) (2.2) 

The equations of homentropic fluid motion are most easily written 
in terms of the enthalpy and velocity: 

+ dlv ^ = 0 (2.3) 

cl 

•V 

- -grad h + vV 2 v (2.4) 

where we retain the viscous term in the momentum equation (2.4) 
and assume that the kinematic coefficient of viscosity is a con- 
stant. In general, the variations in all background or mean 
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thermodynamic and transport properties are of order M^ and will 
be neglected in much of the following work. For the moment we only 
consider the viscosity to be constant. 

B. A KINEMATIC DEFINITION OF SOUND 

To develop a rational theory of aeroacoustics, it is essential 
to adopt a definition of "sound." Goldstein (Ref. 7) has made the 
point that many of the arguments that have raged in the modern de- 
velopment of the subject over what is the "source” of sound are 
rather pointless because there is no common definition of "sound." 

A preoccupation with the question, "what is the source of sound?", 
has in fact led to a host of "exact" but "incomplete" theories that 
are only applicable to the first order problem of aerodynamic sound 
production. For example, consider the Lighthill formulation 
(Ref. 2) where the density is regarded as the primary acoustic 
variable. Lighthill takes four equations (continuity and momentum) 
with five unknowns (pressure, density and velocity) and combines 
them into one equation with five unknowns. The left-hand side is 
the classical wave equation for the density while the right-hand 
side contains all five unknowns and is assumed to be the "source"; 
i.e., the Lighthill stress tensor. While the equation is "exact," 
it is a single "incomplete" equation for five unknowns and only 
becomes useful after several additional hypotheses are introduced. 

The most important is that "sound" is a by-product of the fluid 
flow, the "Lighthill hypothesis." The Lighthill stress tensor is 
assumed to be "known" either by calculation or measurement of tur- 
bulence properties in the flow that are by hypothesis independent 
of the sound that it produces. 

While the Lighthill formulation and hypothesis have been the 
most practical means of calculating flow noise, it excludes by 
decree any consideration of the Interaction of sound with the 
primary flow. A complete theory of sound must be able to cope with 
three fundamental questions: 

1) How is sound produced by a "primary" fluid flow? 

2) How is sound processed by a flow? 

3) How does sound affect the "primary" flow? 

The last question is specifically excluded if one adopts the Light- 
hill hypothesis. 

Many definitions of sound can be stated that are consistent 
with classical acoustics. For example, the density fluctuation 
(or saturation) is a convenient acoustic variable and was used by 
Lighthill. The experimentalist might say that "sound" is that part 
of the pressure fluctuation that is propagating with the local speed 
of sound. Two point time delayed measurements are required to identify 
sound" and it is virtually Impossible to implement such a definition 
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in a theory. We know that the pressure fluctuation is not a valid 
definition of sound because hydrodynamic flows have large pressure 
fluctuations that are nonacoustic. Such considerations led 
Blokhintzev (Ref. 8) to make a distinction between sound and 
pseudosound. Later, Ribner (Ref. 3) split the pressure into two 
parts and identified one as the pseudosound or unsteady hydrodynamic 
pressure that would develop in the absence of compressibility. 

The remaining part is "acoustic" and satisfies a wave equation or 
what Ribner termed a dilatation equation. An important physical 
observation that Ribner stressed was that for sound to be pro- 
duced by or propagated in a flow, there must be local fluid dilata- 
tions . 

Prom classical acoustics we know that sound is a wave motion 
in which the kinetic energy is associated with the local motion of 
fluid particles while the potential energy is stored by local 
isentropic compression of fluid particles. Prom the continuity 
equation (2.3), we observe that a local fluctuation in enthalpy must 
be balanced by a volume change; i.e., 

div V = -ij (2.5) 

a 

2 

If the isentropic compressibility (1/a ) is zero, the enthalpy 
variation acts as a Lagrange multiplier that is only required to 
balance the local acceleration in the momentum equation, but it 
cannot lead to volume changes (dilatations) or sound. Herein, lies 
the difficulty with adopting a thermodynamic quantity as the primary 
acoustic variable. The question always remains as to what part of 
the local fluid motion Is "acoustic." To circumvent this con- 
ceptual difficulty, we start out with a purely kinematic definition 
of sound. 

Suppose we are given the solution of (2.3) and (2.4) subject 
to some initial and boundary conditions. We ask, how can we 
identify that part of the velocity field that Is "acoustic"? From 
(2.5) we first note that there must be local volume changes for 
any sound to exist. Furthermore, these volume changes must be un- 
steady. We introduce an acoustic (or dilatation) potential 4> via 
the Poisson equation 

v2 H ■ It dlv * < 2 - 6) 

and define the acoustic particle velocity 

v a = grad <f> (2.7) 

Given the right-hand side of (2.6) by measurement or calculation, 
we can identify a portion of the velocity field that is associated 
with unsteady volume changes. By our definition, sound is synon- 
ymous with "unsteady compressible flow. " If the notion of a time- 
averaged mean flow Is a meaningful concept we can replace (2.6) by 
the equation 
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2 

= drv (v - <v>) 

(2.8) 

Let u 
Then 

denote the remaining nonacoustic part of the 

velocity field 



V = U + V 

a 

(2.9) 

and we 

can easily prove the 

following relations : 




<v> = <u> 

(2.10) 


curl 

v = curl u - 03 

(2.11) 



div u = 0 

(2.12) 


All of the mean flow (compressible or incompressible) is part of 
y. and all^of the vortical flow (steady and unsteady) is part of 
u . The u field is the entire flowfield less the unsteady com- 
pressible flow. 


It is interesting to note that in the modern theory of turbu- 
lent flow (Ref. 9), it is usually assumed that the turbulent 
velocity field is^incompressible even though the mean flow is com- 
pressible. The u field we have introduced is consistent with 
this concept of a turbulent flow. Our definition of <j> also 
agrees with the usual notion of sound in a homogeneous acoustic 
medium where v a is the acoustic particle velocity and u = 0 . 
The acoustic potential satisfies the classical wave equation 


1 3 <J> 

2 -.2 ~ 
a at 


V c <j) = 0 


and the perturbation acoustic pressure is 


(2.13) 


P’ 



(2.14) 


where p and a are constant . 

The acoustic particle velocity v a ^_s typically very small 
compared to the "primary" flow velocity u . For example, a 140 DB 
sound wave has an acoustic particle velocity of approximately 
.6 m/sec. Typical turbulence velocities in a flow capable of pro- 
ducing a 140 DB noise level would be of the order of a hundred 
meters per second. We use this observation in the following 

development to linearize with respect to the acoustic potential 
while we retain the complete nonlinear primary flow. We remark 
that the linearization is not essential and in fact a nonlinear 
wave equation for <j> can easily be derived that would permit one 
to study the onset of wave steepening and the formation of shock 
waves . 
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We point out that the acoustic potential <j> that we use 
throughout the present work differs from the velocity potential 
that we used in Ref. 1 in an important way. The velocity potential 
was used as the primary acoustic variable in Ref. 1, even though 
it contained the potential flow associated with vortex motion. In 
the solution of specific problems, it was found necessary to admit 
discontinuities in <p . The acoustic potential that we use here is 
physically more appealing and is a smooth function of space and 
time . 


Our acoustic potential is related to the definition of sound 
introduced by Howe (Ref. 10). For a homentropic irrotational 
medium, Howe observes that 3<f>/3t is a measure of variation in the 
total enthalpy since the quantity 

+ H = (2.15) 


is a constant (the Bernoulli constant). Beyond this point, Howe’s 
theory makes no further use of the potential 4> nor does he 
separate the velocity field. Rather he chooses the total enthalpy 
(denoted by B ) as the primary acoustic variable and proceeds to 
derive another ’'exact" but "incomplete" acoustic equation. 


Finally we remark that the need for separating the velocity 
field into acoustic and nonacoustic parts was recognized by Morfey 
(Ref. 11) who investigates the very difficult issue of the acoustic 
energy balance. To distinguish the acoustic energy from the total 
energy in a moving rotational flow, one must identify the acoustic 
velocity and pressure or enthalpy. An important energy theorem is 
discussed below in Section II D. 


C. ACOUSTIC THEORY 


We now develop equations for <}> and u . First, we separate 
the acceleration into acoustic and nonacoustic parts ; i.e.. 


W - if + eraa f- - * X « 


at 

= a| + grad ( 2f + 1 S ra § ♦ ■ £ ) - grad * x 


0) 


where 


(2.16) 


D _ 9 
Dt ~ 8t 


+ U 


grad 


(2.17) 


is the substantive derivative following the nonacoustic motion of 
the fluid. The first term in Eq. (2.16) Is the nonacoustic ac- 
celeration and Is the only term in the complete absence of sound. 
The second term is the particle acceleration associated with 
acoustic motions. The last term in Eq. (2.16) is a Coriolis 
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acceleration that is the essential coupling between the acoustic 
and vortical modes in a free flow. The importance of this term 
will become clear in the subsequent development. 

We substitute Eq. (2.16) into Eq. (2.4) and obtain 

= -grad f-f + vV 2 u + grad cp x u> (2.18) 

where 

H = h + || + .lg. a ^ *1 - - vV 2 <t> (2.19) 


Next, we solve (2.19) for h and substitute the result into Eq. 
(2.3), to obtain a wave equation for <f> : 


1 5L 

T Dt 
a 


( 


D± 

Dt 


grad <j) 


- v V‘ 


DH 


7? Dt 
a 


T=-r- + div 


u 



But div u is not a function of time (see (2.12)). 


( 2 . 20 ) 


In fact. 


div u 



= - 0(M 2 |grad | ) (2.21) 

3 , 

so that the right-hand side of Eq. (2.20) can be replaced by the 
unsteady part of — . Next, we assume (see discussion on p. 10 ) 

that a 2 Dt 


| grad <p\ « |u| 


( 2 . 22 ) 


and linearize Eq. (2.20) with respect to <f> . Also we omit terms 
of 0(M 2 ) on the right-hand side of (2.20). The final set of 
equations of our homentropic acoustic theory can be summarized as 
follows : 

Compressible Primary Flow 


_ 2 , 

- ■ 

i dm , 1 dH\ 

(2.23) 

- V <j> 

T Trt - \ m / 


a 

3 3 


Du 

Dt 

-grad H + v V 2 u 

+ grad <p x w 

(2.24) 

< 


u = 0 

(2.25) 
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with the pressure given by 


fip = -P ||- + P SH (2.26) 

If the primary flow u is incompressible in the mean, the time 
average terms can be eliminated in the foregoing system to obtain: 


Incompressible Primary Flow 


_ V 2 d> - -Hr ~ V' 

? Pt 2 * ? Dt 


1 DH 

~7l Dt 
a 


>+■ _ 
~r = -grad // + vV 2 ui 

div u = 0 


+ grad <f> x to 


(2.27) 


(2.28) 

(2.29) 


The pressure is given by Eq. (2.26) and the speed of sound is a 
constant. We shall be concerned primarily with the last set of 
equations in the present report. 

Our decomposition of the velocity field has led to a natural 
decomposition of the thermodynamic state. The pressure consists 
of two parts; i.e., the acoustic pressure 

6p a = “ p §t (2.30a) 

and the hydrodynamic pressure 

6p = pS H (2. 30b) 

Prom (2.28) and (2.29), we observe that H is the enthalpy asso- 
ciated with the incompressible flowfield and \ \ > H satisfy the 
equations of a hydrodynamic flow (terms in the dashed box of ( 2 . 28 ) 
and (2.29)) with the exception of the Coriolis acceleration term, 
grad <J> x to , in the momentum equation. Thus, in a homentropic flow, 
H is the pseudosound pressure Introduced by Blokinchev (Ref. 8) 
and used extensively by Ribner (Ref. 3). The quantity H is 
dimensionally an enthalpy and we have previously referred to H as 
the "Bernoulli Enthalpy." It is important to note that the H field 
introduced in Ref. 1 has a different definition than in the present 
work. Therein, we used H to generalize the Bernoulli constant of 
homentropic Irrotational flow. The Bernoulli enthalpy is absolu- 
tely constant, except where vorticity is present. Like the acoustic 
variable Introduced in Ref. 1, H is usually discontinuous in a 
specific problem although the physical pressure given by Eq. (2.26) 
is continuous. We use the term "Bernoulli enthalpy" for H in 
the present report with this important distinction. 
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Our acoustic theory is "complete” although the convective 
wave equation for <J> is not exact; i.e., we have linearized the 
acoustic field - a nonessential approximation. The equations for 
the primary flow ( 2 . 28 ) and (2.29), are exact. We can pose any 
of the three fundamental questions of aeroacoustics (see p. 8 ). 

There is no mysterious "source" of sound in the context of the 
present theory. We have by definition isolated the "acoustic" 
and "vortical" (or hydrodynamic) modes (see discussion in Ref. 6 ), 
and have obtained a set of equations that illustrate the essential 
coupling between the two. Energy can flow from one mode to the 
other. We regard the subject of aeroacoustics as the study of 
these mode interaction processes. Whenever we use the term "source" 
in the following discussion, we mean that energy is being converted 
from one form to another. 


D. ENERGY THEOREM 


We prove an energy theorem that brings out the essential 
coupling between the basic modes. Consider the total kinetic 
energy of the primary flow in a fixed volume V , i.e.. 


/ 


u 


Then 


E = I g- dV 


V 


(2.31) 


dE 

dt 


"J S * 


3U HV 

w av 


V 


/ 


dV 


V 


- v 


j^div^-^ u^+ v grad u 2 / 2 ^ 

- grad * • (u x S )1 ( 2 . 32 ) 

3x J 3x° J 


or using the divergence theorem. 


dE 

dt 


= //[-V" + « 2 /2) + v h u2/2 ] 


dS 


- v 


/' 


V 


Wi 3ui 

k 8x J " 8x^; 


IdV 


-/ 


grad <J> • (u x m)dV 


(2.33) 


V 
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We assume that V Is sufficiently large that the surface integrals 
in (2.33) vanish. Thus, we obtain 

( lu 9u \ dv _ f gpad 4, . (u x w)dV (2 . 3*0 

\ 3x J 3x J / J 

V V 

The total kinetic energy of the primary flow can be changed by 
two mechanisms : 

1) Viscous dissipation 

2) Acoustic-vortex interaction 

The conversion of kinetic energy into heat by viscous dissipation 
is an important physical process, but ife not of direct concern to 
the present discussion. In fact, we cannot rigorously account for 
this energy since we have assumed a homentropic model of the flow. 
The main point of our theorem is to illustrate the importance of 
the Coriolis coupling term in the momentum equation and the con- 
version of energy from the primary flow to sound or vice versa. If 
the expression 

grad * (u x w) (2. 35) 

is identically zero throughout the volume V , then no acoustic 
energy conversion can take place. Energy conversion can only 
occur in those regions of the flow where there is a Coriolis ac- 
celeration S x iJ that is nonzero and .where the vorticity w is 
nonzero. It is the work done by the Coriolis force against the 
acoustic particle velocity that results in a local transfer of 
energy between modes. If the energy conversion is primarily in one 
direction, we can "loosely" interpret this as a production process 
where one mode is the "source" of the other. An important corollary 
to our basic theorem is that mode coupling cannot occur in regions 
where the primary flow has a potential. For example, energy trans- 
fer can only occur in the core regions of a vortex flow. 

The energy theorem we have presented Is of quite a different 
type than the results of Morfey (Ref. 11), who is concerned with 
consistent definitions of acoustic energy and fluxes that can be 
used to calculate noise. Our result may not be of direct computa- 
tional value but serves to Illustrate a necessary condition for the 
local conversion of energy. We remark, however, that the acoustic 
formulation we have presented is a good starting point for the 
derivation of energy theorems, following the approach of Morfey. 

E. THE KINETIC "ORIGIN" OF SOUND - THE LIEPMANN ANALOGY 

The foregoing energy theorem has focused on the importance of 
vorticity for mode coupling. Without vorticity and in the absence 
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of finite boundaries, there can be no hydrodynamic flow. It is the 
production of vorticity by the action of viscosity at a boundary 
that is ultimately responsible for the unsteadiness of the primary 
flow and therefore of any aerodynamic sound that is produced. For 
the moment we shall be concerned with the question, "given a 
vorticity field at some Instant of time, what is the mechanism 
whereby sound Is produced in the subsequent development of the 
flow?" To simplify the discussion, we neglect viscosity and con- 
vection in (2.27) and Invoke the Lighthill hypothesis. Then, 


1 3 2 4> „2, _ 1 DH 



-grad H + 


W 2 u 


(2.36) 


div u - 0 


(2.37) 


The hydrodynamic flow can be calculated independent of the acoustic 
field (by the Lighthill hypothesis) and we have the nonlinear 
vorticity transport equation 



grad 


u 


T7 2+ 
vV 10 


( 2 . 38 ) 


and the Poisson equation for the Bernoulli enthalpy. 


V 




8 2 u 1 u J> 

3x i 3x^ 


(2.39) 


From a given initial state, the solution of (2.38) will 
evolve Independent of the acoustics. The Bernoulli enthalpy can 
be calculated by integration of the momentum equation or the 
Poisson equation (2.39). This enthalpy field is the body force 
distribution or acceleration potential required to maintain con- 
tinuity of the Incompressible flow. It acts like a Lagrange 
multiplier and plays no dynamic role In the evolution of the flow. 
The effect of a very slight compressibility is first realized in 
(2.36); that is. 


P 


— (-) 
3h V 0 ) 



(2.40) 


is the Isentropic compressibility of the medium where 1/p is the 
specific volume. The fluctuating hydrodynamic enthalpy produces a 
local unsteady volume change that Is proportional to 1/a^ . It is 
the change of the acceleration potential or Bernoulli enthalpy 
following the hydrodynamic fluid element that is locally responsible 
for the local conversion of hydrodynamic or turbulent kinetic 
energy into sound. The extent to which the local "dilatation" 
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sources cancel to give an overall far field sound pattern must be 
found by solution of the wave equation (2.36). 

To complete our description of the kinetic origin of sound, 
we present a physical analogy that was suggested by Liepmann 
(Ref. 12). Consider the pendulum shown in Fig. 2.1 



Fig. 2.1. The Liepmann Pendulum - Physical Analogy of Sound 
Production. 


If the stiffness k of the rod Is infinite we can calculate the 
primary motion of the pendulum with given Initial conditions. This 
primary motion is analogous to the hydrodynamic flow of our present 
formulation. The tension in the rod is proportional to the cen- 
tripetal acceleration but plays no dynamic role in the evolution 
of the motion. It is a Lagrange multiplier analogous to our Ber- 
noulli enthalpy. If the rod Is slightly elastic, the changing 
centripetal acceleration will excite an oscillation in the rod that 
is analogous to the sound field of our acoustic theory. Energy is 
converted from the rotational motion of the pendulum to the elastic 
oscillation of the rod. The coupling is via the Coriolis acceler- 
ation. If we also Include a small amount of damping in the rod, 
we can actually damp the motion of the pendulum by analogy with 
acoustic radiation energy loss. 

The Liepmann analogy Is one of many analogies that one can 
construct to illustrate a very basic physical principle. When a 
slightly elastic rotating body undergoes a nonuniform acceleration, 
there Is always the possibility of converting energy of the primary 
motion Into an internal vibrational mode. The conversion of 
kinetic energy of a hydrodynamic vortical flow into sound is a 
beautiful illustration of this principle. 
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F. HOW IS SOUND PROCESSED BY A FLUID FLOW? 


The second basic question of aeroacoustics that we have posed 
concerns the interaction of an externally applied sound field with 
a primary flow. For example, if we shine sound on a vortex wake 
or Jet, we know from observation that a scattered sound field is 
produced. In Section III, we consider in detail the interaction 
of sound with a primary vortex flow. Here we briefly discuss the 
problem in a more general context to bring out some of the essen- 
tial features of our acoustic theory. 


Consider the familiar problem of sound interaction with a 
two-dimensional mean parallel shear flow (see Fig. 2.2). The 
perturbation hydrodynamic flow can be represented by a stream 
function ip so that 


and 


Ju ( V ) + I _ i 

1U 0 W 1 3y J dx 


= -£u'(y) 


- ItV 1 ^ 


(2.41) 

(2.42) 


where the prime on u denotes ordinary differentiation with respect 
to y . 0 


Scattered 



Fig. 2.2. Acoustic Interaction with a Mean Parallel Shear Flow 
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Rm + 


(2.46) 


or the single third order equation of Lilley 


- V 2 h 


+ 2u- if*ii 

o 8x9y 


In (2.45) and (2.46), w is the vertical component of velocity 
(see Pig. 2.2). If the viscous terms in (2.43) are omitted, one 
can recombine the results to obtain the Lilley equation (2.47). 


The acoustic and perturbation vortical modes are very clearly 
separated in (2.43). We can discuss the interactions, that result 
when sound impinges on the basic shear flow. The last equation for 
the stream function (2.43) is the Orr-Sommerfeld equation with a 
right-hand side that is the curl of the Coriolis coupling (see 
(2.28)). For most velocity profiles u 0 ('y) the Orr-Sommerfeld 
equation has unstable modes and we expect that the acoustic input 
on the right-hand side will excite these modes. In other words, 
sound can catalyze the process whereby mean kinetic energy is con- 
verted into turbulence. The enthalpy associated with the additional 
turbulence will produce more sound through the first two equations. 
Of course, we cannot use the linear equations to study this 
mechanism in detail because of the inherent exponential growth of 
the unstable modes. 
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Instabilities of the linear perturbation equations have 
precipitated some discussion over the validity of the Lilley equa- 
tion. In practice, any unstable modes that may exist are simply 
suppressed in the construction of the Green r s function of the 
Lilley operator (Refs. 13 » 1*0 so that the matter is primarily 
one of philosophy. If the unstable hydrodynamic modes are going to 
be suppressed (an ad hoc procedure), we can argue that an alterna- 
tive procedure would be to suppress the perturbation hydrodynamic 
mode altogether in (2.43); i.e., set ip = 0 and solve for <j> and 
H (see Section III below). Prom this author 1 s point of view, 
a more serious objection can be raised with regard to the Lilley 
model and its application to the study of acoustic interactions 
with a fully developed turbulent shear flow. There does not appear 
to be any logical way of distinguishing turbulent fluctuations in 
the primary flow from the coherent hydrodynamic fluctuations in- 
duced by sound except when the wavelengths and frequencies of the 
two are significantly different. In practice, this means that 
scattering theories for turbulent shear flows based on the Lilley 
equations are probably only valid in the compact limit. 

G. HOW DOES SOUND AFFECT A FLUID FLOW? 

Our third basic question, the converse of the Lighthill 
hypothesis, is probably the most difficult to examine theoretically. 
The reason is that the uncertainty in calculating a meaningful 
primary flow is usually greater than the effect of any incident 
sound field. Yet, the question has been of considerable interest 
to the experimentalist for many years. The stability of jets to 
external disturbances was discussed by Lord Rayleigh (Ref. 15). 

More recently. Brown (Ref. 16), Hammitt (Ref. 17), Crow (Ref. 18) , 
Bechert and Pfizenmaier (Ref. 19), and others have systematically 
studied the behavior of jet flows when subjected to broad band 
noise and pure tones. Pure tones can amplify broadband noise, 
and the jet itself can be excited by its own sound field or 
externally applied sound. These observations indicate that sound 
is not always a weak by-product of a primary flow as the Lighthill 
hypothesis asserts. Significant feedback mechanisms are possible 
and should be isolated experimentally and studied theoretically. 


In Section V, we show how a simple vortex flow can in fact be 
excited by sound. The basic acoustic feedback mechanism is the 
Coriolis acceleration term in the momentum equation (2.28). While 
this term may be extremely small, it can have a profound long- 
term effect on the primary flow. 

H. COMPARISON WITH OTHER THEORIES 

To conclude our theoretical development, we compare our basic 
equations with several of the acoustic theories that have been 
proposed since the original work of Lighthill. A direct comparison 
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with a single acoustic equation with ’’source" cannot be made for 
reasons discussed earlier. We have a "complete" interactive theory 
without any source. To provide a common basis for discussion, we 
adopt the Lighthill hypothesis and create a "source" model (or 
acoustic analogy) as we did in Section II E. Our acoustic equa- 
tion is (2.36). Our "source" is the hydrodynamic pressure or 
Bernoulli enthalpy that satisfies the Poisson equation (2.39). 
Alternatively, it can be calculated from the solution of (2.37). 

Our equations are formally analogous to the dilatation model of 
Ribner (Ref. 3), although Ribner uses an equation for the acoustic 
pressure instead of the dilatation potential. Also, for low 
speed flows, Ribner would replace the convective derivative of 
H in (2.36) by simple time derivatives; i.e.. 


1 

a 



71 


1 9/V 

2 3t 

cl 


(2.48) 


V 2 w= - a2uluJ 


3x 1 3 x J 


(2.49) 


with the pressure given by 

P 




(2.50) 


sound pseudosound 

If we differentiate (2.48) with respect to time and use (2.49) and 
(2.50), we recover Lighthill ’s equation for p’ ; i.e.. 


3 2 p’ 
a 2 3t 2 


V P ’ = P 


3 2 u i u J ‘ 
° 3x 1 3x J 


(2.51) 


The formal equivalence of our equations to those of Lighthill and 
Ribner in the low speed limit is thus proved. 

An important point must be made concerning the approximation 
of replacing the substantive derivative in ( 2 . 36 ) by the partial 
time derivative in (2.48). There is a sharp distinction between 
the substantive derivative of H and the substantive derivative of 
$ . In an unsteady hydrodynamic flow, both terms in the convective 
derivative contribute equally to DH/Dt . A good example of this 
is the spinning vortex pair considered in Section IV. Only at 
very low speeds (vortex Mach number less than 0.1) where the 
acoustic wavelength Is many, many times the vortex separation 
can one consider the local derivative of H to be a suitable 
approximation of VH/Dt In the acoustic calculation. On the other 
hand, refractive effects in the acoustic operator can be neglected 
for relatively much larger Mach numbers. The importance of this 
point will become more clear when we consider the detailed calcul- 
ation of the sound produced by a vortex pair. 
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The development of Powell (Ref. 20) on vortex noise is re- 
lated to the present work. Where we have called attention to the 
Coriolis coupling between the acoustic field and the primary flow 
Powell has developed a source theory starting with the Lighthill 
formulation that focuses on the Coriolis acceleration of the 
primary flow as the important source term. We can illustrate the 
main points of Powell T s theory with the Lighthill equation (2.51). 
The right-hand side is rewritten as follows 


9 2 u 1 u j ‘ 

8x i ax J * 


div 


D u 
D t 


div(L) + v 2 u 2 /2 


(2.52) 


-> -> 
L - -u X 03 


where 

The solution of (2.51) is of the form 


t = P ° r [div(L + grad u 2 /2)] v >T -> 

" W J |S-y| 


(2.53) 


where the integral is over all space and the numerator is evaluated 
at the retarded time 


t = t - ]x-y|/a . 

Integrate by parts in (2.53) and evaluate the far field to get 


p ' (x) = 


Wa 


/ 


3t" ■ ^ ^ 


Po 8 2 


/■ 


+ ^ |[u 2 /2]dy (2.510 

4-rra 2 x 3t 2 1 
o 1 1 


where 


L x 


i\ 


/ ’-y x 

(u X to) 


(2.55) 


is the projection of the Coriolis acceleration in the direction 
from the source to the observer. For very low Mach numbers, Powell 
argues that the second term is one higher order in Mach number 
compared to the first term. An argument used by Hardin (Ref. 21) 
in an alternate development of the Powell theory is that if the 
region of the primary flow is compact then differences in retarded 
times can be neglected and 
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J [u 2 /2]dy = (E)* (1 + 0 (M) ) (2.56) 

where E is the total energy of the primary a§d the asterisk 
denotes evaluation at the retarded time t - |x|/a . But, to 
lowest order E is a conserved quantity (see Section II. D) so 
that the time derivative in (2.5*0 eliminates the term. Follow- 
ing Hardin (Ref. 21) we expand the first term in (2.54) in a Taylor 
series around the common retarded time and note that the total 
dipole moment (force) is zero in a free flow. The final form of the 
Powell result is 


P’ 



x 1 x J ' 



I 


y^jdy 


(2.57) 


The second time derivative of the moment of the local Coriolis 
acceleration is the dominant source of the far field in the compact 
limit. The Powell theory is particularly easy to implement when 
concentrated line vortices are the dominant source. (See work of 
Hardin, Ref. 22. The main point is that Powell has derived a 
"compact" theory. If phase variations over the "source" region 
(the primary flow) are significant then the arguments leading to 
the simple result (2.57) are not valid. In our discussion of the 
vortex pair (a problem originally considered by Powell) we show 
that phase variations do become important for surprisingly low 
Mach number. In such cases it appears that one cannot avoid a more 
detailed integration over the complete hydrodynamic primary flow. 


We have discussed in previous sections the connection of our 
theory with Lilley (Ref. 4) and with Howe (Ref. 10). The foregoing 
discussion of Powell's theory also pertains to the theory of Howe 
in the form he uses to discuss noise production. In fact, he 
refers to the source term in his (compact) acoustic equation as 
the Powell dipole and writes. 


^ - vV = 

at 


-P Q div(u x to) 


C2.58) 


The importance of compactness in the use of Lighthill's theory 
was discussed by Crow (Ref. 23). In modern applications of the 
theory, as developed say by Ribner (Ref. 3), or Mani and Balsa 
(Ref. 13, 14), it is customary to assume locally compact sources 
and include phase differences between the various localized source 
regions. For turbulence generated noise- at moderate speeds, this 
is probably correct. If the turbulence Mach number exceeds 0.1 
the local source regions themselves tend to generate less noise 
than the compact M/ law indicates, (see Section IV). These non- 
compact effects should be included in high speed noise theory. 
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III. INTERACTION OP SOUND WITH STEADY VORTEX PLOWS 

A. PLANE WAVE - SCATTERING BY A VORTEX WITH CORE STRUCTURE 

We consider the scattering of a plane wave propagating along 
the x-axis from a single vortex with a radial core vorticity dis- 
tribution ft(r) , as shown in Pig. 3.1 



Figure 3.1 Scattering of Plane Waves by a Vortex 


The tangential velocity field is given by 


V(r) 



(3.D 


We assume that the maximum Mach number of the vortex is sufficien^ 
tly small that we can neglect quadratic terms. Our basic inter- 
action model becomes 




- V"d> = 


1 _ 9 H 

2 at 


(3.2) 


V 2 H = div (grad <f> x kft) (3.3) 

Since H is of order (M) we can neglect the convection terms in 
Eq. (3.2). Now let 


2k 



where 


-iwt 

4) = Re(<f>^ + <f> s )e 




ikx 


(3.3) 

(3.4) 


To lowest order (the Born or first scattering approximation) the 
equations for the amplitude of the scattered field are 


where 


♦b + = 

- S 1 

(3.5) 

II 

£ 

CM 

t> 

div 

(3.6) 

2ikV 9c ^i 

ik .. 


ar 9 0 

- IT" 

(3.7) 

= (grad 4> 1 


(3.8) 

•H 
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•H 

t*r-3 
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II 
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The first term in S-. is due primarily to the interaction of the 
sound with the potential flow "mantle" of the vortex while S 2 is 
the direct Coriolis interaction with the vortex core. 


While we can in principle, solve for <j> and H everywhere 
in the flowfield, we are primarily interested in the scattered far 
field. Thus, we solve Eq. (3.5) in the far field to obtain 


i 

~ IT 


H^ ) 1) (kr) 


/■ 


-ikr. 


-> 

■y 


s 1 (y)dy 


(3.9) 


where 


The solution for 



the H field is 


H{ y) 


fel 

1 ln|y - z 

J 

11 

f(y - z) 

2ir 

1 ■+■ **>.2 

J 

|y - z ! 


| div ^ 2 dz 

•S 2 dS 


( 3 . 10 ) 


(3-11) 


where we have used the divergence theorem to obtain the last result. 
We calculate the scattered fields due to the mantle and core as 
follows. Let 


= 4> c + 4> m 


(3.12) 
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where 


, _ k . rr( 

♦o " 55: H o 


‘’(wJV” 1 


-ikr,. -y . _ 

Hi y)ay 


>m Hsr H o 1) ( te) 


/ 


-11? -J, 


I % 

7 se / 


(3.13) 

(3.14) 


Substitute Eq. (3. 11) into Eq. (3.13) and reverse the order of 
integration to obtain 


k „(: 
Bua n o 


_ . C -ikr ^ f -ikr-, «y £ ^ 

1 ^(kr) I e 1 3 2 dz* Je ^ dy (3.15) 


The last integral is easily evaluated. We get 




H^Ocr)?^ j”( 


-+• -> 
-ikr 1 « z 


S 2 aS 


(3.16) 


Now substitute Eq. (3.8) and Eq. (3.4) into Eq. (3»l6) and use polar 
coordinates to obtain 


k<j>. ( 1 ) 

H^ ±; (kr) sin 0 | e 




-iky[cos (0 - B T ) - cos 0’] 

fl( y)dy d0' 
(3.17) 


The integration over 0 f can be carried out with the substitution 

0» = e" + 0/2 

so that «» 

6 = - ~ 4>?H ^ 11 ^ (kr) ) yftJ ( 2ky sin 0/2)dy*sin 0 (3.18) 

T c 2a 1 o j o 

To evaluate cj) we introduce polar coordinates in Eq. (3.14) 
and substitute for ^i from Eq. (3*4). Integrate once by parts 
with respect to 0 f and the integral becomes 

-iky[cos (0 - 0 T ) - cos 0'] 
e Vdy d0 r (3.19) 


k_ 

2a 


<{kH^ (kr) f* 
m 2 a Y i o 9e J 

The integral over 9’ is the same as In Eq. (3.17). Thus, we get 

2 

k 77 j.o„(l) I --WT (2ky sin 0/2)dy cos 0/2 (3.20) 


m 


(kr) 

a i o 




Integrate by parts with respect to y to obtain 


<j> m = fa ^i^o"^^ 1 ^ j*ynJ Q ( 2 ky sin 0/2 )dy cot 0/2 (3.21) 
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where we have used Eq. (3.1) in the differential form 


fy yV - y£2 (3.22) 

By comparing Eq. (3.20) and Eq. (3.21) we note that the mantle 
scattering has been reduced to an Integral over the core vorticity. 
This Is not too surprising since the mantle velocity field Is de- 
termined uniquely by the core vorticity. The Interpretation as 
mantle scattering rather than core scattering is important however. 
Observe that Eq. (3-21) is singular in the forward scattering 
direction indicating a strong focus. This is due to the great 
extent of the flowfield around a single vortex that presents a very 
large cross-section to plane waves. Because of this focus effect 
we cannot evaluate the total sound power. Also, the basic scat- 
tering theory must be regarded as invalid in the forward scattering 
direction. To get a better understanding of this focus effect we 
consider a line sound source in the next section. 


We combine Eq . (3.21) and Eq. (3.18) to obtain the final result 
for the scattered field. 


_ rk 
~ TFa 


(kr) cos 9 cot 


r s c (e) 


(3.23) 


where 


and 


s o ( 0 ) 
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yftJ o (2ky sin 0/2)dy 
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(3.24) 


(3.25) 


The intensity of the scattered field normalized to the incident 
intensity is 


where 



<2 

<p 



3t 


at 




S = cos 0 cot 


0 

7 


s c ( 0 ) 


(3.26) 


(3.27) 
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Concentrated Vortex 


For a concentrated point vortex the core factor S c is unity 
and the angular distribution of the scattered intensity is 

S 2 = cos 2 8 cot 2 | (3.28) 

For a core without mantle the same factor is 

S 2 = sin 2 6 (3.29) 

a result previously obtained by Howe (Ref. 10), and Yates and 
Sandri (Ref. 5). The basic directivity pattern, Eq. (3.24), agrees 
with a previous result of Muller and Matschadt (Ref. 24), and is 
compared qualitatively with the core alone in Fig. 3.2 



Figure 3*2. Qualitative Comparison of Basic Directivity 

with Core Scattering Only 

The dipole character of the core scattering is in sharp contrast 
to the quadrupole pattern with strong focus in the forward direc- 
tion when the mantle is included. Further discussion of the core 
results is given below. 

The scattered intensity is proportional to the factor 
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where X is the wavelength of the incident sound and 


(3.30) 


is the characteristic acoustic length scale of the vortex. 

\ << X the vortex is essentially transparent to the incident 

sound. For X Y = 0(X) or greater the validity of the Born approx- 
imation may be T in question unless the core compactness is such that 
the amplitude of the scattered field remains bounded. We consider 
now the effect of various core structures. 


(3.3D 

If 


Solid Body Rotation 

If the vorticity is constant inside the core radius r c then 


= 


TTT, 


0 < r < r 


= 0 


r > r 


(3-32) 


and we readily calculate 


S = 
c 


J^(2kr sin 0/2) 

kr sin 0/2 
c 


= 1 - 


(kr c ) 


(3-33) 


sin 9/2 , kr c << 1 


Exponential Core 

For an exponential distribution of vorticity 
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2-rrrf 
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-r/r e 

e 


(3.34) 


we calculate 
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P 2 -3/2 

S c - [1 + (2k r 0 ) sin 0/2] 


= 1 - 6(kr ) 2 sin 2 0/2 , kr << 1 (3-35) 


Betz Core 


Pop aircraft wing vortices, it has been shown from the theory 
of Betz (Ref. 25), that the vorticity is distributed according to 
the formula 


ft =3-£p|- 
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The core scattering factor is 
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(3.36) 
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sin 0 / 2 j dt 
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^ - 2) sln 0/ " 2 


kb 
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For any core we note that the principle effect is to reduce the 
scattered intensity. Scattering in the extreme forward direction 
is less affected except for large values of the compactness para- 
meters (kr ) . 


To compare the relative importance of the three cores, we 
require that the total vortex strength and polar moment of vorti- 
city be the same for each distribution. Then the three core lengths 
are in the ratio 



= (1.5 , 0.5 , 4.65) 


(3.38) 


The reduction in scattered intensity due to the core is given 
approximately by the first term in the asymptotic form of each 
formula for S c . These terms are in the following proportion 


(Solid Body, Exponential, Betz) 


9 


(1.125 



1.5 , 2 . 16 ) (3.39) 
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The solid core is the most effective scatterer while the Betz core 
is the least effective. In general, if the vorticity is more dis- 
tributed, the amplitude of the scattered sound field will be less 
than for the equivalent concentrated vortex. 

The angular distributions of the scattered intensity for the 
three core types are compared in Pigs. 3.3 a,b,c. For the expon- 
ential and Betz cores the effect is a uniform reduction of the 
scattered intensity. For the solid body core we note a more compli- 
cated diffraction pattern as the core compactness ratio is increased. 


Solid Body Core Without Mantle 

We point out that the previous calculations of Howe (Ref. 10), 
and Yates and Sandri (Ref. 5), for the scattering from a finite 
solid body vortex core without mantle are in error. In both cases, 
only the scattering due to the Coriolis interaction with the core 
is considered. Direct refraction in the wave operator was neglected. 
We can use the present results to obtain the correct answer. The 
vorticity distribution Is 





r) 


( 3 . 40 ) 


where r is the total vortex strength of the core neglecting th.e 
shell of vorticity at r = r Q . The scattering factor for this 
core is 


0 

S = cos 0 cot 2 J£(2kr sin 0/2) 

(kr ) 2 

= — — sin 20 , kr c << 1 (3.41) 

The pattern is that of a pure quadrupole in contrast to the dipole 
pattern reported earlier (see Fig. 3.2). Also the efficiency of 
the core as a scatterer is less by the compactness factor (kr ) 2 /8. 
Refraction due to direct interaction with the velocity field in 
the core is just as important as the Coriolis interaction. 

B. ACOUSTIC INTERACTION WITH DISCRETE WEAKLY INTERACTING VORTICES 

The results of the previous section are useful for estimating 
the importance of core structure. The basic scattering pattern is 
singular, however, and it is difficult to use the results in a 
straightforward calculation of the scattered Intensity. In the 
present section, we consider localized noninteracting (or weakly 
interacting) vortices in the plane and a more general sound field. 
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OLID BODY CORE 
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ody core on intensity distribution of scattered plane waves. 




EXPONENTIAL CORE 
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Effect 
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Later we consider the particular case of a line source from which 
we obtain the general scattering Green’s function. The approach 
in this section is somewhat different in that we adopt a Lagrange 
formulation of the basic interaction process. 

Consider a single line vortex interacting with an arbitrary 
acoustic field in the plane (Fig. 3.4). The velocity field is the 
sum of two parts 



Figure 3.4 Interaction of an Arbitrary Sound Field with a Line 

Vortex 


= u(x - t.) + 


v a = grad <f) 

k x y 

u(y) = y 2 = grad x 


(3.42) 


where t. is the Lagrange coordinate of the vortex, i.e., 


g|=v a (A, t) (3.43) 

and 

V = I_ 

T 27 r (3.44) 
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where 

terms 


T is the vortex strength. The H field is calculated in 
of the velocity potential x > i.e.. 


H = 





(3.45) 


Again, we consider only first-order interactions and neglect all 
terms quadratic in y . Then 
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Our acoustic equation becomes 
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(3.47) 


The term on the right-hand side is the dilatation rate produced by 
the acoustic acceleration of the vortex. The coupling term on the 
left-hand side is again due to convective refraction by the poten- 
tial flow mantle. Both terms are of the same order of importance 
to the interaction process. 


We remark, without proof, that the scattered sound field cal- 
culated with Eq. (3-47) and incident plane waves is identical to 
the result obtained in the previous section (see Eq. (3-23)), with 
the core factor set equal to unity. The calculation is straight- 
forward and was carried out as a check on the Lagrange formulation. 


The extension of Eq. (3.47) to N noninteracting vortices is 
immediate. We simply add up the direct interactions on the right- 
hand side and use the convective velocity field due to all of the 
vortices on the left-hand side. We get 
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\ Mr + % £, 

a 3t a n=l 


- *).!-* - V 2 , 


3x9 1 
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a n=l 


n 1 3x3t/ J 


(3.48) 


n 


36 



r 


where A Is the Lagrange coordinate of the nth vortex, the notion 
of nonin^eract ing vortices can be made more precise with Eq. (3*^8). 
The ’s are considered to be slowly varying functions of time 

by comparison with the time required for a sound wave to traverse 
the complete collection of vortices. For example, if we think of 
the aircraft vortex wake, individual vortex filaments do not move 
relative to each other when viewed in a coordinate plane moving 
with the aircraft. 


C. LINE SOURCE INTERACTION WITH A LINE VORTEX 


An important special case of the foregoing theory is that of 
a line source of sound interacting with a line vortex as shown in 
Fig. 3 . 5 . For simplicity we neglect any back reaction of the vortex 
flow on the acoustic source. In application, the effective source 
impedence would have to be estimated or derived from a separate 
calculation. 




Line monopole source Vortex at origin 


Figure 3.5 Line Source Interaction with a Line Vortex 

The results can be used to synthesize a more general scattering 
Green’s function for multiple vortices. Also, we can use the 
results directly to estimate the effect of vortex wakes on engine 
noise . 

For a simple harmonic source at x = x Q we have 

_ m \ -imt 

+i = (lei* - \l )e (3.49) 

and the equation for the scattered amplitude Is 
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r k x x 

u _ 2n i+i 2 


(3.50) 


(3.51) 


We use the general results of Section IXI-A, Eqs. (3*9 ~ 3*14) to 
write 


» =0 + (}> 
s Y c m 


(3.52) 


where 
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From Eq. (3.49) 


^ k - |! H< x) (ki? - x i ) 

3y 1 ly - x Q | 


(3-53) 


(3.54) 


(3.55) 


and 


P 1 ) = *5 

■»y 4> |x ol 


(3.56) 


Consider the core scattering. Substitute Eq. (3.51) and Eq. (3.56) 
into Eq. (3.53) to get 


■c “ fe H o 1} ( kr) 


r-ikV? r_l4l.*o^o H U) (tato)d 5 

•/ y n 


= - (j>° (kr ) H| 1) (kx Q ) sin 0 (3-57) 

This is the familiar dipole pattern that we obtained with plane 
waves (see Eq. (3.18)). The vortex core sees a local plane wave 
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whose amplitude is reduced by the distance factor (fkx ) . We 

note that the second Hankel function in Eq. (3.53) is the°complete 
function while the first function is to be replaced by its far 
field asymptotic form, (see Eq. (3*10). 


The ’ calculation of the scattered sound due to the mantle is a 
slightly mor^ complicated problem than the core calculation. We 
use Eq. (3.5*0 and introduce polar coordinates. Then 


-+■ 

u* 




inn - 


(3.58) 


and after Integration by parts in 0* we get 

kr 3 l,' lky cosCe - 0,) H ( 1 ) 
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(3.59) 


To carry out the integration we expand the Hankel function in the 
Integrand in a partial wave series (see (Ref. 26), page 827). 


where 


H o 1)(k ly - 5 ol> =Y2 e (-i) m G m (y >x o ) 

m=-°° 

G m ( y» x o ) = J m (ky)H m 1)(kx o ) y < x o 

" J m (kx o )H m 1)(ky > y > x o 


( 3 . 60 ) 


(3.61) 


With formulas in Ref. 27, page 254, the integral in Eq. (3.59) can 
be evaluated in detail. The final series solution for the total 
scattering is summarized below 
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Y J 
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(3.65) 


where the argument of all of the Bessel functions in the definitions 
of T , A , B is kx For m greater than kx we remark 

that Whe s®ries m converge§ rapidly and is easy to calcfilate. 


The important parameter that enters the point source calcula- 
tions is kx Q . It is 2ir times the number of wavelengths between 
the acoustic source and the vortex center. In general, this para- 
meter is not small as we shall see in specific application to vortex 
wakes. 


The foregoing results are for point vortices. We can correct 
for finite core structure with the results obtained in the previous 
section. We simply multiply Eq. (3.62) by the core factor defined 
by Eq. (3.24). Also, we write down the result for multiple vortex 
scattering : 
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(3.66) 

(3.67) 


( 3 . 68 ) 


where ft (y) is the vorticity distribution of the nth vortex core 
and T n n is its total strength. Also 
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(3.69) 


The core correction we have introduced is a good approximation pro- 
vided the source center is several core radii from the vortex center. 
The foregoing results can obviously be used as a Green’s function to 
calculate the scattering of more complex source distributions. 


Far Field Intensity 

We calculate the far field intensity for the single vortex- 
single source combination. We have 


where 
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ImtJ> = (k| x - x Q |) + H^ 1) (kr)C(0) 


-iwt 


and 


r 
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(3.71) 


(3-72) 


We use the vortex center as the reference origin (see Fig. 3.5) and 
expand the source in the far field 


n 'i -v hi ikx^ cos e 

Hq 1 (k| x - x Q | ) = H^ 1) (kr)e 

The final result for the far field intensity is 


(3.73) 
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where 
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?[< 


(ImC) sin (kx Q cos 8) 


+ (Re C) cos (kx Q cos 0 




I 1 “ P a 
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k (<f> ± ) 

irkr 


(3.74) 

(3.75) 


is the source intensity. We compare the overall intensity with the 
source intensity in DB via the formula 


A db = 10 log 10 (3.76) 

1 

Note that A DB can be positive or negative. Typical results are 
plotted in Figs. 3.6 a,b, and c for the source location kx D = 10 
and three values of the amplitude factor Tk/a . Note that the back 
scatter is virtually zero for the smaller values of T k/a . Max- 
imum scattering is in the forward direction at approximately 30 de- 
grees off the x-axis. The scattering pattern is reflected in the . x- 
axis, if the sign of V is reversed. The strong forward scattering 
may be compared with the singular focus in the case of plane waves. 

D. SCATTERING OF ENGINE NOISE BY AN AIRCRAFT VORTEX WAKE 

We now apply our vortex scattering theory to estimate the 
effect of an aircraft vortex wake on engine noise. We conclude our 
investigation with a specific application to the DC-9 in a typical 
take-off configuration. Strong flap vortices with r of 
133 mvsec are located approximately 5.4 m from the engine. The 
engine frequency at the peak Strouhal number of 0.3 is about 240 Hz. 
The acoustic wavelength at this frequency is about 1.5 m. The essen- 
tial dimensionless parameters that we need are 

kx o = 25 

and 


We assume that the left engine is scattered by the left wing vortex 
and the right engine is scattered by the right wing vortex. Also, 
we assume that the engine sources are incoherent so that we can 
superimpose the intensity in the far field. The interactive noise 
field for this particular configuration is shown in Fig. 3-7. We 
estimate that the engine noise in the aircraft ground signature is 
enhanced by 3 or 4 DB due to the close proximity of the flap vor- 
tices. For higher frequencies in the engine noise spectrum there 
could possibly be greater amplification although the vortex core 
structure would have to be accounted for. 
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d Intensity. 



Pig. 3.6c. Directivity of Intensity Relative to Unscattered Intensity. 
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3 DB 3 DB 


Fip;. 3-7- Enfine flap-vortex interactive noise for DC-9 in 
take-off configuration. 
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These simple estimates suggest that significant variations of 
ground engine noise signatures can be caused by the trailing vortex 
configuration. If measurements are to be made of engine noise 
ground signatures it is recommended that a standard flap configura- 
tion be specified. These interference noise estimates could also 
explain discrepancies between measurements made with the same engines 
on different aircraft. Finally, we remark that detailed calculations 
of specific engine and wake configurations can be carried out with 
the multiple vortex scattering theory we have presented. It would 
be desirable, however, to develop the counterpart of Eqs. (3.66 - 
3.69) for a point source near a line vortex. 

In conclusion, we remark that we have focused on the weak 
interaction problem in the present investigation, i.e., the inter- 
action parameter Tk/a must be of order one or less. It Is equally 
important to consider large values of this basic parameter which 
means large wave numbers or high frequency. Ray acoustics theory 
would be the logical theory to apply. In this regard we remark 
that in recent work of Burnham et al (Ref. 28), results of ray theory 
have been used to design a wake vortex tracking system. The device 
works at a nominal frequency of 3*5 kHz so that for typical aircraft 
vortices (T = 500 mVsec) , the interaction parameter is of the order 
of 100 . The backscatter is weak but the device works quite well at 
short range. It is currently operational at Kennedy airport where 
it is used to detect vortex decay during normal airport operation. 
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IV. PRODUCTION OF SOUND BY VORTEX FLOWS 


A. THE COROTATING VORTEX PAIR 

We first consider an elementary two-dimensional vortex flow to 
which we can apply our basic acoustic formulas. The problem is that 
of two vortices of equal strength that spin about their centroid as 



The problem has been investigated in various ways by different 
authors (Refs. 20,29) and is an excellent problem for comparing 
various aeroacoustic theories. 

Source Evaluation 


We use the production formulation of equations ( 2 . 36 ) and (2.37) 
with viscous effects omitted. Thus 


1 

2 

a 



V cf> 


1 DH 

2 D t 
a 


Du 
D t 


- grad H 


div u = 0 


(4.1) 

(4.2) 

(4.3) 
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The hydrodynamic velocity field that satisfies (4.2) and (4.3) is 
most easily expressed in terms of a velocity potential. Thus 

u = grad x 



(4.4) 


with 
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y 
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The enthalpy 1-1 
formula 


is most easily calculated from the Bernoulli 
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The final result is 
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where 



a = 2(9 - <j>) (4.8) 


The quadrupole character of the H field is evident from the 
appearance of the double angle in (4.7). For later reference we 
note the asymptotic behavior of H ; i.e. 3 
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(4.9) 



The first term is the steady state pressure due to the effective 
single vortex of strength 2 T at the origin. The remaining terms 
are quadrupole, octupole and higher order components of the unsteady 
hydrodynamic pressure field. 


The acoustic "source" with the present production model is 
the substantive rate of change of the H field. After a tedious 
differentiation we obtain the following result: 


where 


DH _ 

Dt 


— P(R,a) 
r 


o 


(4.10) 


F(R,a) = R 2 (R 2 -D(R 2 -7)sln a (4.11) 

( 1+R -2R cos a) d 

It is important to point out that in the calculation of H and 
DH/Dt that the convective derivatives are as important as the local 
time derivative. Only in the far field is 3#/3t a valid approx- 
imation of DH/Dt . 


The form of the source function F(R,a) is plotted in Pig. 

4.2. The axes are node lines along with the circles with normalized 
radii R = 1 and /f . There are three distinct annular quadrupole 
regions that are alternately 90 degrees out of phase - the far 
field R > /J s the region 1 < R < /7 and R < 1 . We turn now 
to the calculations of the noise generated by each of these regions 
and the interesting phase cancellation between them. 


The Acoustic Problem 


The general solution of the acoustic equation (4.1) is given 
by (Ref. 26, p. 893). 
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We Introduce the normalized variables 


a ■ x/r 0 

3 0 = y/r 0 (4.14) 
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and the eddy Mach number 


M = 



ftr 
o 

a 


(4.15) 


where S2 is the angular velocity of the spinning pair. Then 
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The integrals in (-4.16) can be evaluated if we note that F(R,a) 
(see (4.11)) is a periodic function of a and has a Fourier series 

i . e . , 


F(R,a) = 1 S (R) sin na 


where 
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The coefficients S (R) can be evaluated explicitly with known 
results of Ref. 30, p. 113. We get 

S (R) = n — ^ ( - 7 ^ -I , 0 < R < 1 

1 + R 


= n(R -7) 
R 2n (l+R 2 ) 


R > 1 


(4.19) 


Substitute (4.17) into (4.16) and carry out the integration over 
t to obtain 
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We remark that the last expression could be used to evaluate $ 
and all acoustic properties in the near and far field. The results 
could then be used to calculate refractive effects and higher order 
acoustic-fluid interactions. In the following, we evaluate the 
acoustic far field only. 

For R >> 1 we have 


(2Mn|^-^ o | ) * (2MnR)e“ 2;LMnRoCOS ^ e ^ (4.21) 

where it Is understood that the Hankel function on the right-hand 
side Is to be replaced by its asymptotic expression (see 3.10). 
Substitute (4.21) Into (4.20) and carry out the integration over 
<J> . We get 


$($,T) ^ Re l (-l) n H (1) (2MnR)e 21n( ‘ G MT) C (4.22) 

n=l 2a*V ° n 

o 


where 


n 


00 

/ 


C„ = I S v ,(R)J 2n (2MnR)RdR 


n 


(4.23) 


The total sound power radiated is given by the expression 
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0 n=l n 

The far field solution is a series of outgoing partial waves, 
the leading term of which is a quadrupole whose frequency Is twice 
the basic rotational frequency of the vortex pair. The amplitude 
C n of each partial wave must be evaluated numerically. The weight 
of the integrand in (4.23) gives us a measure of the virtual source 
region of the flow (see Fig. 4.2). For very small Mach number and 
n not too large the weight of the Integrand Is In the far field 
where R = 0(l/2Mn) . We replace S n by its asymptotic value and 
obtain in the compact limit 
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The first three values of 


are as follows : 
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n_ 
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2M 2 /3 
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(4.26) 


The strength of each higher order multipole decreases with an 
additional M 2 factor 

From (4.24) and (4.25) we observe that in the compact limit 

P ^ 16tt^ ( pa^r Q )M^ (4.27) 

7 

The sound power radiated varies as M' in accordance with known 
results for„compact two-dimensional aeroacoustic theory (see Ref. 31). 
A simple power law was obtained for the spinning vortex pair 

by Powell (Ref. 20) because he considered a segment of a three- 
dimensional ring pair. The M? law was obtained by Muller and 
Obermeir (Ref. 29) who solved the problem in the compact limit 
by matched asymptotic expansions (see below). 

In the compact limit all acoustic theories seem to converge 
to the same answer for the acoustic power radiated. We are now in 
a position to place bounds on the validity of the compact limit 
for this simple problem. The multipole coefficients (4.23) and 
the sound power were evaluated numerically over a range of Mach 
numbers. The results are presented in Fig. 4.3. The M? law is 
given by the dashed line and the numerically calculated total sound 
power is given by the solid line. The various multipole contribu- 
tions are also shown separately for comparison purposes. The total 
sound power is pure quadrupole for eddy Mach numbers, less than 
0.3 • The higher order poles contribute to the total sound power 
for higher Mach numbers although the basic calculations we have 
carried out are suspect for these Mach numbers. An interesting 
mathematical point is that the multipole series seems to diverge 
for M > 2/e where e is the base of the natural logarithms. 
Numerical calculations also indicate that this is the case. 

The real significance of our calculation is for M < 0.3 . 

First, we observe that the compact limit is asymptotically valid 
for Mach numbers less than 0.1 . As the Mach number increases 
toward 0.3 the power radiated diminishes by 15 DB due to the non- 
compactness of the vortex pair structure. We remark that the 
acoustic wave length is 30 times the radius of the vortex pair 
at M = 0.1 ! This result gives an idea of how large the wave 
length must be to treat an eddy as compact . In a high speed tur- 
bulent flow we would expect eddy Mach numbers of the order of 0.2 
to 0.3 • Our results indicate a significant reduction in the sound 
power due to eddy noncompactness. 
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Matching in the Compact Limit 


A useful and familiar technique in aeroacoustics is that of 
matched asymptotic expansions. Crow (Ref. 23) used the technique 
to make a general critique of the Lighthill theory. We remarked 
earlier that Muller and Obermeir (Ref. 29) used the method to solve 
the spinning vortex pair problem. It is a simple and powerful 
tool for calculating acoustic fields in the compact limit. We 
conclude our discussion of the vortex pair by using the idea of 
matching and illustrate the limitation of the method to the 
lowest order compact limit. 


The asymptotic expansion of the hydrodynamic pressure field 
of the vortex pair is given by (4.9). We know that this field 
is transformed into an acoustic pressure field for sufficiently 
large R . In some intermediate overlap region the far field 
approximation of fi should match the near field approximation of 
the acoustic field. We carry out the matching for the leading 
term in each multipole that we obtain from (4.9); i e. s let 
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The far field acoustic pressure (enthalpy) with the same phase as 
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where A is to be determined by matching. We expand (4.29) for 
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Now equate (4.30) and (4.28) and solve for A 
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We get 
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Substitute (4.31) into (4.29) and expand for large R to obtain 
the acoustic far field. 
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where 
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2n-l M 2n-2 


2 ( 2n-l )~ ! 


(4.33) 


The last result is Identical to (4.22) with C n given by the 
compact approximation (4.25). The matching procedure yields the 
correct results for the multipole coefficients in the compact limit. 
We have not been able to obtain corrections for noncompactness 
via the matching procedure or any other method. It appears that a 
detailed integration over the near field pressure must be carried 
out to obtain these corrections. 


B. JET IMPINGEMENT NOISE 

In recent experimental work of Preisser and Block (Ref. 32) the 
noise of subsonic jets impinging at normal incidence on a plane 
wall has been carefully measured. In earlier work of Snedeker and 
Donaldson (Ref. 33) extensive measurements of mean flow and tur- 
bulence properties of impinging jets were carried out. The acoustic 
data indicate noise levels 10 to 15 DB greater than the noise of 
the free jet. The flow data Indicate that the turbulence levels 
near the impingement point are not substantially different from 
those in the free jet. What is the mechanism for the excessive 
noise? By simple imaging of the free jet sources in the plane 
boundary, we can argue that the noise should be 6 DB greater than 
the free jet. What is the origin of the other 5 to 10 DB? 

From our discussion of the general problem of noise production 
in Section II, the physical mechanism for the impingement noise 
is evident. Consider a turbulent eddy (for example, the two-vortex 
model) that is convected along a streamline near the impingement 
point. The acceleration of the eddy due to streamline curvature 
is a real noise production mechanism. The local ^/-field will be 
magnified by the acceleration of the mean flow and so DH/Dt will 
be greater. To estimate the order of the magnification, we con- 
sider the Poisson equation for H ; i.e., 

V 2 H = -(uV) ± . 

3 3 O 

= -2u 1 .u J 'l - u 1 lu J '! (4.34) 

tJ J 1 3 3 S 1 

where we now use tensor notation and introduce a mean u^ and fluc- 
tuating velocity field u^ . Recall in C2.49) and (2.51) that the 
source of H is the same as the source of Ribner f s pseudosound 
equation and Llghthill's acoustic equation. The magnitude of the 
source of the ^/-field is a direct measure of the magnitude of the 
sound field. We estimate the increased impingement noise by 
evaluating the source term, in particular the shear noise, in 
(4.34) for a free and impinging jet. 
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The geometry of the impingement model is shown in Pig* 4.4. 

We suppose that the impingement point is in the fully developed 
turbulent region of the free jet. The mean impingement flow field 
can be well-approximated with a general result of Barnes and Sullivan 
(Ref. 34) who considered a Gaussian velocity profile impinging on an 
infinite plane wall. Unfortunately, their exact solution is a series 
of hypergeometric functions that converges very slowly except in 
the vicinity of the impingement point. There it can be shown that 
the mean flow can be approximated by a stream function that in 
cylindrical coordinates is of the form 

= Cr 2 z (4.35) 

The constant C is a measure of the curvature of the mean flow as 
we shall see below. 


The velocity components are readily evaluated with (4.35) 
and we get 
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(4.38) 


are the nonzero ChristoffeL symbols for cylindrical coordinates. We 
substitute (4.36) and (4.38) into (4.37) to obtain the following: 
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u = C 


= -2C 


(4.39) 


where all other components of u . vanish. Note that the sum of 
the three terms in (4.39) is zeri^in accordance with the incom- 
pressible continuity equation. Also, we have introduced the phys- 
ical components of the mean velocity in place of tensor components. 
If we substitute (4.39) into (4.34) it can easily be shown that 
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Fig. 4 . J| . Geometry of Jet Impingement Problem. 
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(4.40) 


V 2 ^/- 6C ■»— + (u 1 u^* ) . . Impinging 

t f ’ aJ Jet 

Shear Self 

noise noise 

Following RIbner, we call the two terms shear noise and self noise. 
We point out that the basic physical mechanism of the shear noise 
term is the "acceleration" of turbulent eddies in the curved stag- 
nation point flow. The resulting acoustic radiation is very much 
analogous to the phenomenon of Brehmstrahlung radiation in atomic 
physics as pointed out by Williams (Ref. 35). 

To bring out the essential role of the impingement acceleration 
we compare (4.40) with the corresponding result for a free jet 
where we model the mean flow as a parallel shear flow. We readily 
obtain 

1 2 H = -2 1^- -f- + (u 1 V ’ ) , , Free (4.4l) 
3r 3z ,i,j Jet 

The experimental results of Snedeker and Donaldson (Ref. 33) in- 
dicate that the magnitude of the turbulence is not appreciably al- 
tered by the impingement process. The turbulence is essentially 
frozen along streamlines until several radii from the impingement 
point where it re-accomodates to the wall shear. On the other 
hand, the experiments of Preisser and Block (Ref. 32) indicate that 
the acoustic source is localized within a few radii of the impinge- 
ment point. Thus, we conclude that the enhancement noise must 
result from the shear noise term. In (4.40) the coefficient C is 
given by ( see (4.39)) 

c = -§£ (4.42) 

We compare C with 3w/3r in (4.4l). It seems plausible that 
3u/3r_ for the impinging jet should be the same order of magnitude 
as 3w/3r of the free jet. The length scale for both terms is 
the radius or diameter of the free jet. The velocity change in 
both terms must scale with the average mean flow velocity. Assuming 
that the two terms are the same order of magnitude, we observe that 
the ' impingement source is greater than the free jet source by a 
factor of 3 • Without carrying out the detailed acoustic integrals, 
we can see that the factor of 3 Increase In the source strength will 
translate into a factor of 9 in the sound power or approximately 
10 DB. If this is added to the 6 DB due to imaging of the free jet, 
we have an estimated noise increase of 16 DB in rough agreement 
with the acoustic over pressure measurements of Preisser and Block. 

C. A SUGGESTED PROBLEM 

A detailed calculation of the radiated noise should be carried 
out with the source term (4.40) and measured or calculated jet 



turbulence data. Time did not permit this calculation under the 
present contract. While It is important to carry out the foregoing 
calculation it is equally important to illuminate the physical 
mechanism that leads to the increased noise. A two-dimensional 
model problem that would illustrate the basic enhancement mechanism 
is the following: We have carried out a detailed calculation of 

the noise generated by two co rotating vortices in a free flow. 

The next step is to imbed the pair in various accelerating mean 
flowfields; e.g., the stagnation point flow or mean shear flow. 

The noise field should be enhanced by the mean acceleration and 
become directional. If the eddy Mach number is sufficiently low, 
we can use the compact Powell theory (Ref. 20) as extended by 
Hardin (Ref. 21) to calculate the noise. For high speed jet im- 
pingement, we should account for noncompactness of the eddy vortex 
pair as shown in Section IV of this report. For the axial jet, 
the noise due to impinging vortex rings should be calculated. The 
work of Davies t et al. (Ref. 36) would again be the easiest to 
apply in the compact limit. 


V. EXCITATION OF A FLUID FLOW BY SOUND 
A. DISCUSSION OF EXPERIMENTAL RESULTS ' 

It has long been recognized experimentally that sound can 
have an effect on the development of a fluid flow. For example, 
the transition from laminar to turbulent flow in a boundary layer 
can be altered at will with the incident sound field (Ref. 37). 

This is not surprising since the basic laminar flow is known to be 
sensitive to small disturbances of any kind for sufficiently large 
Reynolds numbers. 

More interesting is the fact that turbulent flows can be 
stimulated or excited by impressed sound fields or by their own 
sound. For example. Brown (Ref. 16) has observed the excitation 
of jets by pure tones over a range of frequencies. He found, for 
example, that eddies form in the outer part of the jet at the 
frequency of the incident sound. For certain excitation frequencies, 
there was almost an explosive excitation and spreading of the jet. 

In other experiments, Hammitt (Ref. 17) has recorded the noise of 
a jet and played it back to the same jet. The result was a large 
amplification of the turbulence in the jet. Crow and Champagne CRef. 
18) considered the excitation of jet turbulence by upstream dis- 
turbance pure tones and carefully measured the enhancement of the 
turbulence at the excitation frequency and its harmonies. More 
recently, Beckert and Pfizenmaler (Ref. 19) have shown that pure 
tone can amplify the jet broadband noise. All of these results 
indicate that a free turbulent flow can be modified significantly 
by sound. Because of these facts, there has been a renewed 
interest in the excited jet as a diagnostic tool and by those 
who are interested In noise suppression (see Ref. 3.8 e.g.). 
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If the turbulent structure of a jet could be controlled acoustical- 
ly, there is the possibility of designing a more efficient noise 
suppressor . 

B. THE LIEPMANN ANALOGY REVISITED 


In Section II E, we used the Liepmann analogy to illustrate 
the basic mechanism of noise production by an unsteady hydrodynamic 
flow. With a slight extension of the analogy we can Illustrate the 
principle of fluid excitation by sound. Suppose that we excite 
the vibrational mode of the Liepmann pendulum by some external 
source. The vibrational energy can then flow into the rotational 
motion of the rod. By analogy, if we shine sound on a rotating 
fluid element, we should be able to excite the fluid motion. In 
the subsequent discussion, we carry out a detailed calculation of 
the spinning vortex pair and show that this is indeed the case. 

The following example illustrates the direct transfer of 
acoustic energy into the hydrodynamic vortical mode. An equally 
Important mechanism to be considered is the catalytic effect of the 
sound field on the production of turbulence from the mean flow. 

We mentioned this point earlier in reference to acoustic stimula- 
tion of boundary layer transition. Also, we have seen from the 
linearized equations (Section II F) how the Coriolis coupling 
can lead to acoustic excitation of the basic instabilities of a 
laminar flow and, therefore, of the transition process. The analysis 
of transition phenomena is difficult without the complication of an 
acoustic field and is beyond the scope of the present study. We 
simply call attention to their importance and consider a more simple 
model problem that contains part of the basic physics. 

C. EXCITATION OF THE COROTATING VORTEX PAIR 


We consider two equal corotating vortices with a plane wave 
incident sound field as shown in Fig. 5.1. The Lagrangian equa- 
tions for the vortex centers t and B are 


at 

dt 


icx(ft-S) 

, ->■ , ? 

| A-B | 


_a. -r -r 

+ iv sin(ki*A-wt ) 

Cl 


(5.1) 


d§ 

dt 


r*- . . 

kx(B-A) 

I A-B | 


+ iv sin(ki-B-o)t) 

cl 


(5.2) 


where v is the magnitude of the acoustic velocity. We introduce 
coordinates for the center of vorticity and the relative spacing. 


a = M 


r = J-g 


Equation (5*1) and Eq. (5.2) become 


(5.3) 
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(5. ft) 


= ?v a sin(ic? ) cos 

§ = 2y + 2?v & cos(^iUt)sin C5.5) 

The acoustically induced motion of the center of vorticity leads to 
scattering. We discussed this problem in detail in Section III. 

In the following discussion, we fix the center of vorticity and 
consider the acoustic effect on the relative motion. The error 
incurred is second order In the amplitude of the acoustic field. 


We Introduce the following scaling and dimensionless notation: 


2r‘ 


r = 2r o (X,Y) 


= 


2 “ r o 


kr. 


t = 


y M J 

The equations for X and Y become 


Y 

k = kr 


(5.6) 


dX 

dx 


+ e sin <X cos fix 


R 


dY _ X_ 
dT r 2 


where 


^ o o 

= X + Y 


(5.7) 


(5.8) 



plane wave 



Pig. 5.1. Acoustic Excitation of Vortex Pair. 
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and 


e 


2r v 
o a 


Y 



(5.9) 


is the ratio of the acoustic velocity to the orbit velocity of the 
vortex. An essential feature of the excitation in (5-7) is that 
there must be a phase difference between the velocities seen by 
each vortex. If the wave number k tends to zero, the incident 
sound looks like a hydrodynamic field to the pair and no excita- 
tion can result. 

We now assume e << 1 and consider the perturbation form of 
(5.7). Let 


X = cos x + eX’ 
Y = sin x + eY ’ 

Then to lowest order the equations for X' 


( 5 . 10 ) 


are 


X T - X T sin 2 t + Y’ cos 2 t = sin(K cos t)cos fix 
Y 1 ~X r cos 2x + Y* sin 2x = 0 

Two independent solutions of the homogeneous equations are 


(5.11) 


X^ = cos x + 2x sin x 
Y| = sin x - 2x cos x 


X£ = -sin x 
Y£ = cos x 


(5.12) 


Note that the first solution has a linear secularity; i.e., for 
arbitrary initial conditions the perturbation solution will grow 
linearly in time. With the homogeneous solutions (5.2), we can 
write down the general solutions of (5.11) in the form 

x. 


X f = (cos x + 2x sin x) 
x 


/ 


cos t sin(xr cos t)cos fit dt 


+ sin x 


/ 


(sin t - 2t cos t)sin(< cos t)cos fit dt 
x 


(5.13) 


Y f = (sin x - 2x cos x) 
x 


/ 


cos t sin(K cos t)cos fit dt 


/ 


- cos x I (sin t - 2t cos t)sin(< cos t)cos fit dt (5.1*0 
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where we have assumed zero initial conditions. For small k and 
large x it can be shown that the integrals in (5-13) and (5.14) 
are bounded except when ft is 2, in which case the perturbation solu- 
tion grows quadratically in time; i.e., 

kt 2 

X T 'v. — £ — sin x 

4 2 
V* T 

Y 1 ^ cos T (5.15) 

From the perturbation analysis, we expect that the nonlinear sys- 
tem will have a maximum response at the acoustic frequency comes- 
sponding to ft = 2 . Note that this is also the first acoustic 
frequency of the vortex pair alone (see Section IV). From the 
scale relations (5.6) we also note that at the resonance condition 

k = ka = 2 77 j- = 2M 


or 


A_ _ TT 

a M 


( 5 . 16 ) 


For a Mach number of order 0.1 we can have a tuned resonance con- 
dition with an acoustic wave length some 30 times the radius of the 
vortex pair. The resonance criteria is a frequency matching con- 
dition and does not require that the acoustic wave length match 
the eddy length scale. 


To verify that resonance excitation of the vortex pair does 
occur we solved (5.7) numerically and used the results to compute 
the accelerations and the acoustic source DH/Dt from (4.10) and 
(4.11). We use the latter in the unnormalized form 


DH 
D t 


y 3 r 2 (r 2 -r2)(r 2 -7^)sin(e-.|>) 
r 2 (r^-Sr^r 2 cos 2(e-<(>) + r°) 2 


(5.17) 


where r , <j> are the instantaneous polar coordinates of either 
vortex and r , 0 are polar coordinates of a point in the flow- 
field source region. The principle cause of the acoustic source 
is the acceleration induced by the incident acoustic field. Thus, 
we give a phase plot of X versus Y and a plot of DH/Dt versus 
time for each case. The observation point for DH/Dt is shown 
in Fig. 5.2. We choose M = 0.1 and e = 0.06 for all cases. 

The magnitude of the excitation was chosen somewhat large for the 
purpose of visualization. The important point to observe in the 
following results is the relative magnitude of the response to the 
excitation. For reference, we first give the unexcited acoustic 
source in Fig. 5.2. The acceleration phase plot is a circle (not 
plotted). Note the anharmonic content in the plot of DH/Dt . 

This is the source of the higher harmonics (see (4.26)) in the acoustic 
radiation of the vortex pair. In Figs. 5.3 and 5.4 we give the 
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Unexcited pair 



Pig- 5.2. Acoustic source of an unexcited vort 


Y 



Fig. 5.3. Acceleration Phase Plot for Detuned Acoustic Excitation; 
fi = 1.5 , M = 0.1 j e = 0.06 . 
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acceleration and DH/D t for a highly detuned case ft = 1 . 5 • 

There is very little change from the unexcited case. In Pigs. 

5.5 and 5*6 we present the results for ft = 2. The change in 
DH/Dt is approximately 25 percent of the unexcited amplitude, 
four times the magnitude of the excitation. In Figs. 5-7, 5*8, 
and Pigs.- 5.9 s 5.10, we present the results for values of ft 
slightly more ft = 2.1 and slightly less ft - 1.9 than the 
linear resonance condition. In the first case, the result Is much 
larger values of the acceleration and DH/Dt . Note that (DH/Dt ) max 
is about 50 percent greater than the unexcited case, indicating an d 
amplification factor of about 8 . For ft = 1.9 the magnitude of 
the acceleration and DH/Dt Is attenuated by a factor of about 8 . 
Thus all incident acoustic frequencies In a narrow band around the 
fundamental pitch ft = 2 are capable of amplifying or attenuating 
the basic vortex motion. 

D. QUALITATIVE COMPARISON WITH EXPERIMENT 

The elementary mechanism we have described for excitation of 
fluid motion by sound can be used to explain in part some of the 
experimental observations. In the experiments of Brown (Ref. l6) 
and Crow and Champagne (Ref. 18) pure tones were found to enhance 
turbulence In the frequency band near the excitation. In fact. 

Brown observed the formation of discrete vortices in the edge of 
the jet with an angular frequency equal to the excitation frequency. 
The jet would also become highly unstable for certain excitation 
frequencies. The process of vortex shedding at the nozzle lip is 
also excited by the incident sound and certainly plays a role in 
the subsequent formation of the jet. However, it seems that the 
complete jet can present a large cross section to the incident 
sound. As such, the free flow excitation mechanism we have des- 
cribed can be equally Important . 

We can use the present results to speculate on the mechanism 
of broadband noise amplification by pure tones that was reported 
by Bechert and Pfizenmaier (Ref. 19). Suppose that a jet flow- 
field is composed of "small scale" random turbulent eddies super- 
imposed on "large scale" coherent vortex structure. Further sup- 
pose that the random turbulent motion is acoustically more efficient 
than the large scale structure. Now, If we ensonify the jet with a 
pure tone that is properly tuned, the vortex motion associated with 
the coherent structure can be energized at the expense of the 
acoustic energy in the pure tone. The random turbulent flow is 
subsequently excited by the more energetic vortex flow and would 
be expected to radiate more broadband noise. 
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VI. CONCLUSIONS 

Because of the comprehensive scope of our study, we summarize 
our conclusions for each section separately: 

A. SECTION II - ACOUSTIC THEORY 

Our general theory of acoustics for homentropic flows is self- 
contained and complete. We have shown that the three basic ques- 
tions of aeroacoustics can be investigated systematically with the 
theory. The definition of sound and the concept of Bernoulli 
enthalpy are the fundamental building blocks. The Importance of 
Coriolis acceleration for coupling of the acoustics and vortical 
modes is Illuminated. 

B. SECTION III - INTERACTION OP SOUND WITH VORTEX PLOWS 

1. Por plane wave scattering from a vortex the overall 
directivity is quadrupole-like with a strong focus in 
the forward direction. Our results are in agreement 
with those of Muller and Matschadt (Ref. 24). 

2. A concise formula for core scattering is derived that 
only involves a quadrature over the core vorticity. 

Vortices with distributed core vorticity (like Betz) 
scatter less than a more concentrated vortex. 

3. The overall scatted amplitude is proportional to r k/a . 
Core structure is important If kr c is of order 1 where 
r c is the geometric extent of the core vorticity. 

4. A general theory of acoustic scattering from discrete 
weakly interacting line vortices is presented. Core 
effects are Included using the plane wave result. Por 

a line vortex, the forward scattering focus (obtained with 
plane waves) is eliminated. Scattering is a maximum at 
about 30 degrees from the forward scattering direction. 

The parameter kx where x is the distance from the 
source to the vortex enters the scattering formula as a 
geometric attenuation factor. 

5. The scattering theory is used to estimate the scattering 
of engine noise by wing vortices. For the DC-9, it is 
shown that scattering of Jet engine frequencies near the 
peak Strouhal number can be of order 3 or 4 DB. 

C. SECTION IV - PRODUCTION OF SOUND 

1. An explicit calculation of the noise generated by a 
corotating vortex pair is carried out. For low Mach 
number (M < 0.1) the M' production law is obtained in 
agreement with other two-dimensional compact acoustic 
theories . 
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2. For 0.1<M<0.3 it is shown that the noise is reduced 
by 10 to 15 DB due to noncompactness of the source region 
even though the radiation is basically quadrupole. 

3. It is recommended that the corotating vortex pair be used 
to model jet impingement noise. 

4. An estimate of jet impingement noise is made. We conclude 
that the impingement noise is due primarily to acceleration 
of turbulent eddies in the curved stagnation point flow. 

The "shear noise 1 ' source is amplified by a factor of 3 and 
the sound power is amplified by a factor of 9 or about 

10 DB. The additional 6 DB due to imaging gives a total 
of about 15 DB. These estimates are the order of magnitude 
of Preisser and Block (Ref. 32). 

D. SECTION V - EXCITATION OF A FLUID FLOW BY SOUND 

1. It is shown that the corotating vortex pair can be excited 
or de-exclted by an incident acoustic field. The excita- 
tion is a maximum when the incident acoustic frequency is 
approximately twice the rotational frequency of the pair. 
Amplification factors of 4 to 8 have been calculated. 

2. The simple excitation model can qualitatively explain many 
of the experimental observations with excited jets. 
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